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It was shown recently that, without jeopardizing the success of the A cold dark matter model 
on cosmic scales, the modified Newtonian dynamics (MOND) can be derived as an emergent phe- 
nomenon when axionlike dark matter particles condense into superfluid on the galactic scales. We 
propose in this paper a Dirac-Born-Infeld (DBI) scalar field conformally coupled to the matter com- 
ponents. To maintain the success of MOND phenomenon of dark matter superfluid on the galactic 
scales, the fifth force introduced by the DBI scalar should be screened on the galactic scales. It 
turns out that the screening effect naturally leads to a simple explanation for a longstanding puzzle 
that the MOND critical acceleration coincides with present Hubble scale. This galactic coincidence 
problem is solved, provided that the screened DBI scalar also plays the role of dark energy on the 


cosmic scales. 


I. INTRODUCTION 


Recently, a novel theory of dark matter (DM) super- 
fluidity [1] [2] was proposed to combine the success of 
modified Newtonian dynamics (MOND) [3][5] on galac- 
tic scales with the triumph of the A cold dark matter 
(ACDM) on cosmic scales. The MOND turns out to be 
an emergent phenomenon of DM itself on galactic scales 
due to a MOND-like force between baryons mediated by 
superfluid phonons of the axionlike particles condensed 
as superfluid with a coherence length of order the galac- 
tic size and a critical temperature of order micro-Kelvin. 
The ACDM model is eventually recovered beyond galac- 
tic scales when the fraction of particles in the condensate 
decreases with increasing temperature due to larger ve- 
locity dispersion and hence larger DM temperature in 
galaxy clusters. 

It was known as the galactic coincidence [6] that a 
critical acceleration scale appears in various seemingly 
unrelated Kepler-like laws of galactic dynamics, which 
cannot be simply explained in a common way in the 
context of the cold dark matter (CDM) scenario. How- 
ever, MOND predicts such a universal acceleration scale 
ag © 107 1?^m/s?, which should intriguingly happen to be 
of order the present Hubble scale Hg ~ ag or more boldly 
the cosmological constant scale Af ~ Mg,a2. Although 
MOND now emerges from DM itself on galactic scales in 
the context of DM superfluidity, the galactic coincidence 
still manifests itself as an input parameter in order to 
fix other parameters to their preferred values. It should 
be in any case striking that the dark matter and dark 
energy sectors have such a common scale even though it 
is currently unclear whether it is just a coincidence or 
smoking gun for new physics. 

It was also known as the cosmic coincidence that the 
energy density used to account for the late-time cosmic 
acceleration happens to be the same order of magnitude 
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as the matter components today. Alternative to the stan- 
dard cosmological constant scenario, one might as well 
consider a slowly rolling scalar field known as dynamical 
dark energy (DE) with proper screening mechanisms [7] 
to hide the fifth force from the local tests of gravity. To at 
least alleviate the cosmic coincidence, the energy density 
in the scalar field should at least track [8] [9] the back- 
ground energy density and then grow to dominate the 
energy budget at late times. Either the screening mech- 
anism or tracking behavior can be realized if general in- 
teractions between dark energy and matter components 
are concerned. 


In this paper, we propose a very simple explanation 
for the galactic coincidence problem by conformally cou- 
pling a Dirac-Born-Infeld (DBI) scalar field with local 
matter components. To effectively screen the fifth force 
mediated by the DBI scalar field from the MONDian 
force mediated by DM superfluid phonons on galactic 
scales, the galactic coincidence aj = A?/2gMpi ~ Ho 
is derived, provided that the DBI characteristic scale 
A* ~ MRH ~ (meV)* coincides with current critical 
energy density for conformal coupling g ~ O(1). This 
allows us to interpret the DBI scalar field as a dynamical 
DE in the presence of a conformal coupling term. The 
equation of state (EOS) of our DBI dark energy mimics 
that of Chaplygin gas. 


This paper is organized as follows. In Sec. |II| we review 
the DM superfluidity and define the MOND transition 
scale. In Sec. we propose a DBI-like scalar confor- 
mally coupled with the matter component to solve the 
galactic coincidence problem. In Sec. the possibility 
of our DBI scalar playing the role of DE is explored. The 
final section is devoted to conclusions and discussions. 


II. DARK MATTER SUPERFLUID 


In the nonrelativistic regime, DM superfluid is 
effectively described by the MOND Lagrangian with a 


conformal coupling term to baryons, 
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where DM particle m is of order eV to ensure the forma- 
tion of Bose-Einstein condensation and the phonon exci- 
tation X = 6—m®—(V0)2/2m is described by the Gold- 
stone boson @ for a spontaneously broken global U(1) 
symmetry under the external gravitational potential ®. 
The dimensionless parameter o and dimensionful param- 
eter A can be fixed later by inputting the MOND critical 
acceleration ao in order to reproduce the MONDian pro- 
file. For static spherically symmetric profile 0 = ut-4-(r) 
at constant chemical potential and baryons distribution 
Ty = —py(r), the equation of motion (EOM) 
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can be integrated for the X « 0 branch to obtain 
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for k > u — m with My(r) = 4r Jo r'"2dr! p(T’), which 
admits a MONDian acceleration, 
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if one identifies 
o? A? 
Mp = ao, (5) 


hence a ^ O(1) for A ~ meV. 

The general picture of DM superfluidity is that the 
DM halo core where galaxies are located is almost en- 
tirely condensed and the dynamics is dominated by the 
MONDian force mediated by the DM superfluid phonons, 
whereas galaxy clusters are either in a mixed phase or en- 
tirely in the normal phase just as those on cosmic scales. 
Therefore, it is natural to define a MONDian transition 


radius 
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in the context of the DM superfluid core with core radius 
TMOND containing the total mass of M. To see that this is 
a reasonable definition, consider a DM halo with central 
density po ~ M,,,/rg and core radius ro = \/M,,G/ao; 
one obtains a constant surface density poro ~ M, /rg ~ 
ao/G independent of galaxy luminosity found recently 
by several astrophysical observations [I0][13]. One can 
even reproduce a sort of baryonic Tully-Fisher relation 
(BTFR) [1416] Mr ~ pora ~ (ao/G)rà ~ v*/Gao by 
using poro ~ ag/G and ao ~ v?/rg. The MONDian tran- 
sition radius thus serves as a natural separation between 
the MOND regime r < rg with ay < ag and the Newto- 
nian regime r > rg with an > ag where ay = GM, [r?. 


III. DBIONIC SCREENING 


The action of the scalar field we propose in this paper 
has the form 


SpBIT,, = | atay- ( 9A - A100) 
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which will be referred to as the DBIT;, action for short. 
It should be kept in mind that the same symbol A used 
in our action has nothing to do with that in the 
action (ip, although they actually coincide as we will 
see later. Here, f is the determinant of the Friedmann- 
Robertson-Walker (FRW) metric of a 3-brane moving in 
a five-dimensional Minkowski space with two time dimen- 
sions, 


ds? = —dw? + f,,dz"dz". (8) 


Here, the Gaussian normal transverse coordinate w(x) = 
A~*6(a) is written in terms of the DBI scalar field (x). 
The first term in DBIT,, action can thus be inter- 
preted as a cosmological constant term, 


P =J d'z/-g(-A*) = J dtz,/—f(-A*y~"), (9) 


in terms of the induced metric g,, = fuv — A~*0,,¢0, 
on the brane, and the inverse of the induced metric is 
just gt” = fH” + A~4y20"¢0"o with an abbreviation 
y= 1/VA- ATOY. 

The first term in differs from the standard DBI 
action 
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by a flipped sign in front of the derivative term, which as 
we will see is essential for the so-called DBIonic screening 
mechanism [I7]. It is worth noting that the first term in 
(7) also differs from 
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in standard DBlonic screening by an overall sign of the 
action, which as we will see is also essential for the scalar 
field to mediate a repulsive fifth force and to drive the 
late-time acceleration. The second term in the DBIT,, ac- 
tion describes a conformal coupling of the DBI scalar 
with the trace of the energy-momentum tensor of back- 
ground matter fields with strength g ~ O(1) from the 
stringy perspective. 

Suppose the DBI scalar field ¢(r) with a static and 
spherically symmetric profile is coupled to a static local 
source Tin = —Pm(r); then, the EOM 
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can be integrated to give 
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where a DBI transition radius 
1 gM 1a (14) 
ropi = = | —— 
rn A 4r Mp 


is introduced to separate the DBI regime r > rppr with 
repulsive force 


f = —2g7an (15) 


from the Newtonian regime r < rppy with screened force 


7$ Mp 

To retain the success of the MOND phenomenon of DM 
superfluidity on galactic scales, the DBI force should also 
be screened from the MOND force on the galactic scale, 
which renders an identification of the DBI transition ra- 
dius (14) with the MOND transition radius (6), 
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Therefore, the galactic coincidence 
nu 
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is derived, provided that 
A‘ ~ Mà H2 ~ (meV)! (19) 


for a conformal coupling g of order unity. It turns out 
as a nice surprise that A4 coincides with current critical 
energy density and A in the DBIT;, action (7) matches 
that in the MONDT, action (ip. This is why we use 
the same symbol for the scale A in both actions and 
(7), which shares the same scale with the cosmological 
constant. 


IV. DBI DARK ENERGY 


The repulsive feature of the DBI force and the unex- 
pected match of A* with the current critical energy den- 
sity inspire us to explore the possibility of our DBI scalar 
field playing the role of dark energy. 

We start with the total Lagrangian 


V-FL = V-fLg + V-fLer +v -fLm, — (20) 
where 
£4 = —A*/1— A-4(0¢6)?; (21) 
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Lm = Ln( fur, Y). (23) 


A. Backreaction on matter 


In the absence of the conformal coupling term, the 
matter component is supposed to behave as a pressureless 
fluid with the trace Tm = —pm of the energy-momentum 
tensor T5 = (2//—f)6(/—fLm)/df"”. In the presence 
of the conformal coupling term, the matter field could 
exchange momentum by interacting with the DBI scalar 
field. "Therefore, the conformal coupling term would 
necessarily introduce an effective pressure in the matter 
fluid, and the effective EOS parameter of matter could 
in principle deviate from zero. We will show below that 
such a deviation from pressureless fluid can be made ar- 
bitrarily small for a sub-Planckian DBI scalar. 

'The EOM of the DBI scalar field for a spatial homoge- 
nous profile $(t) is simply 


" T gTm 
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according to the Euler-Lagrange equation 
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In the absence of the conformal coupling term, the 

energy-momentum tensor of the DBI scalar field can be 
computed as 
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with its energy density and pressure of the form 
pa =X; (27) 
pp =— May. (28) 


In the presence of the conformal coupling term, the con- 
servation equation of the above energy-momentum tensor 
should be written as 


gTm 
VES = ^ Mp Ov, (29) 


where the temporal component of the above equation 
reads 


H 9pm 1 
bo + 38H (ps + po) = ux 9, (30) 


which is consistent with the EOM (24). 

In the_absence of the conformal coupling term, the 
EOM (24) has a trivial solution o = 0, and the EOS 
parameter 


ACA (31) 


would simply imply a cosmological constant with wg = 
—1. In the presence of the conformal coupling term, the 


EOM (24) cannot admit such a trivial solution db = 0 un- 
less $ is always equal to zero, which is of less physical in- 
terest. Therefore, our DBI scalar should generally behave 
as a dynamical Chaplygin gas [I8] py = —A°/pg with 
phantomlike EOS parameter and superluminal sound 
speed cà = p/p = y~7, where the closed timelike 
curves are argued to be evaded within the regime of valid- 
ity of the effective field theory (EFT) due to chronology 
protection [20] P1]. With slow-roll condition ¢ < A?, our 
DBI scalar could serve as a candidate for the DE sector. 
We will show below that such a slow-roll condition can 
be satisfied for a sub-Planckian DBI scalar as well. 

To derive the conservation equation for the matter 
component, we start with an alternative definition of the 
energy-momentum tensor for the DBI scalar, 
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with its energy density and pressure of the form 
EV gó : 
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In the presence of the conformal coupling term, the con- 
servation equation of the above energy-momentum tensor 
should be written as 


vererer — 9 9 Ty 35 
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where the temporal component of the above equation 
reads 


. go . 
por + 3H (por + por) = M Pm (36) 


which is also consistent with the EOM (24). 

Since the total energy-momentum tensor is conserved, 
the conservation equation of the energy-momentum ten- 
sor of the matter component is thus 


Tus (37) 


where the temporal component of the above equation 
reads 
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The source term on the right-hand side of above equation 
can be accounted for by recognizing the effective EOS 
parameter of the matter component as 


Pm + 38H pm = — 


Wm = —— - 1. (39) 


Therefore, the backreaction of the DBI field on the mat- 
ter component due to the conformal coupling term can 
be safely neglected in the field region à < Mp; of the 
DBI scalar for conformal coupling of order unity. From 
now on, we will take a fiducial value g = 1 for the con- 
formal coupling in order to solve the galactic coincidence 
problem. 


B. Steady flow assumption 


In the rest of this section, we will work with the as- 
sumption, called the steady flow assumption, that the 
energy flow from the DBI scalar to the matter compo- 
nent is conserved. We define the energy flow as the 
energy-momentum tensor associated with the conformal 
coupling term 


T T . 
T$ = Tir —T$-fwCer: (40) 
then, steady flow assumption is expressed as 


veTéT = S 8 (Ta) — 0, 41 
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where the temporal component of the above equation 
reads 


Pm + dpm = 0. (42) 


The steady flow assumption simply states that, although 
the energy-momentum tensors of the DBI field and mat- 
ter field are not separately conserved as indicated in Eqs. 
and (7). there is no loss during the energy trans- 
fer from the DBI scalar to the matter component and 
the total energy-momentum tensor of the DBI field and 
the matter field is conserved, namely, VEI, TVETD = 
Vere = 0. We will justify numerically the steady 
flow assumption below. 

With the steady flow assumption, one can solve the 
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Mp, go 2go (2) 
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analytically by combining Eq. (38) with Eq. (42), where 
$o = la = ag) with present-day scale factor ag = 1 
and W(z) is the Lambert W function defined by z = 
W(z) exp[W (z)]. Hence, the evolution equation of 
the matter component can be directly integrated to give 
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(44) 
The evolutions of DBI field, the effective EOS param- 
eter of matter component, the matter energy density, 
and the conformal coupling term are presented in Fig. 
[1] The backreaction of the DBI field on the matter com- 
ponent is negligible during the matter dominated era as 
long as a sub-Planckian field value for the DBI field at 
present is specified. However, the effective EOS param- 
eter of the matter component will eventually approach 
—1 in the future, causing an unavoidable vacuum de- 
cay to matter, saving us from big rip singularity as we 
will see. The steady flow assumption is justified by a 
constant conformal coupling term. At small scale fac- 
tor a < 1, the evolution of the Lambert W function 
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FIG. 1. The evolutions of the DBI field, the effective EOS parameter of the matter component, the matter energy density, and 
the conformal coupling term with respect to the scale factor for initial conditions ¢9/Mp; = 107+, 102, 107°. 


W(a?) ~ a? compensates the evolution of the matter 
component pm ~ a-? to render a constant conformal 
coupling term Tm ~ W(a?)pm ~ const. At a large 
scale factor, the constant nature of the conformal cou- 
pling term is nontrivial. 

The evolution of the energy density of the DBI field 
can be solved numerically by rewriting Eq. (30) as 
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With numerical solution p; (a), one can evaluate all other 
quantities like 
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where the effective EOS parameters wF (a) and wor (a) 
of the DBI scalar field are defined by rewriting Eqs. 


and in a form without the interacting term, 


jo + 3H(1 + wg! )pe = 0; (51) 
Por + 3H(1-- wor)por = 0. (52) 


'The evolutions of the above quantities are plotted in 
Fig. The division of DBI fluid from matter fluid 
is somewhat artificial since the DBI scalar and matter 
component are coupled together. However, the differ- 
ence between definitions and of the energy- 
momentum tensor of the DBI scalar are shown to be 
negligible in Fig. therefore, we will just stick to 
Eq. (26) for the sake of simplicity. We also compute 
the evolution of the Hubble parameter by 3MàH (a)? = 
Perla) + px (a) + py(a) and the fractions of energy den- 
sity by Q;(a) = pi(a)/3M3,H (a)? in Fig. [3| It is worth 
noting that the DBI scalar relaxes its phantom nature by 
vacuum decaying to matter, preventing the matter com- 
ponent from being diluted away and leading to a constant 
Hubble parameter in the asymptotic future free of big rip 
singularity. 


C. Slow-roll conditions 


Last but not least, it is the slow-roll condition 
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FIG. 2. The evolutions of the energy density of the DBI field and their effective EOS parameters with respect to the scale 
factor for initial conditions $o/Mpi = 107+, 1077, 107°. 
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FIG. 3. The evolutions of the Hubble parameter and fractions of energy density with respect to the scale factor for initial 
conditions ¢o9/Mp; = 107+, 1077, 107°. 


that allows us to interpret our DBI scalar as a candidate 
for the dark energy sector. To evaluate analytically the 
EOS parameter of our DBI DE, we propose a second 
slow-roll condition, 


? ; <1, € <1, (54) 
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Recalling that the factor y = 1/4/1 + A-49? and the 


matter component Tm = —Pm = —3M?,H?Qm and the 
galactic coincidence A^ = 4g? M2 H2, one can immedi- 
ately derive from the above equation the EOS parameter 
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where the reduced Hubble parameter E = H/Hp is un- 
derstood and the conformal coupling g is surprisingly 
canceled out. Testing Eq. with the present value 
of matter fraction Qmo e 0.3, one finds the present value 
of the EOS of our DBI DE, 
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Woo 1.023, (57) 
perfectly matching the Planck 2015 constraints [22]. A 
distinct feature of our DBI DE is that w¢ọo and Qmo 
are strongly correlated without other free parameters en- 
countered. Although behaving mildly like the phantom 
at present, our DBI DE will relax its phantom nature 
by vacuum decaying to matter, preventing matter from 
being diluted away, resulting in a constant Hubble pa- 
rameter and leading to a de Sitter future free of big rip 
singularity. The validity of the first and second slow-roll 


conditions and is presented in Fig. 


V. CONCLUSIONS AND DISCUSSIONS 


It was recently claimed that the axionlike dark matter 
particles can condense on galactic scales as a superfluid, 
the phonons of which mediate MONDian force between 
baryons, and thus MOND arises as an emergent phe- 
nomenon of dark matter itself. The standard ACDM 
model is recovered on cosmic scales in the presence of 
dark matter particles in the normal phase instead of the 
condensed phase. We have proposed to study the pos- 
sible origin of the MOND critical acceleration scale in 
the context of dark matter superfluidity. We have intro- 
duced a DBI-like scalar field conformally coupled to the 
matter components. It turns out that the MOND criti- 
cal acceleration is roughly at the same magnitude with 
the present Hubble scale, provided that the conformally 
coupled DBI scalar plays the role of dark energy. 
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FIG. 4. The evolutions of the first slow-roll condition o? < 
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Mee with respect to the scale factor for initial conditions 


$o0/Mpi = 1071, 1077, 107°. 


However, one might be concerned with the possible 
ghost problem of our proposal. In canonical quantum 
field theory, a Lagrangian with a wrong-sign kinetic term, 
after canonical quantization, usually admits the nega- 
tive norm states with negative energy, namely, the ghost 
states. If there are no other fields directly coupled to the 
ghost field, it would not cause us any trouble. However, if 
there are other fields with a correct-sign kinetic term di- 
rectly coupled to the ghost field, the vacuum would be un- 
stable because it could generate a pair of ghost particles 
with negative energy and a pair of normal particles with 
positive energy. We argue that the possible ghost prob- 
lem might not be as pronounced as it appears to be due 
to the following three features encountered in our model. 
First, the Hamiltonian density turns out to be positive 
and bounded below, which suggests that there might be 


a stable vacuum where ghost particles can condense. Sec- 
ond, the equation of motion is second order in the time 
derivative, which might evade the ghost problem from 
the view point of Ostrogradsky’s theorem. Third, even if 
the ghosts indeed exist, they are indirectly coupled to the 
matter fields via the trace of the energy-momentum ten- 
sor. Since the matter fields act as a source term, there 
are simply no sources for ghosts to be generated when 
DBI-like scalar field come to dominate. This might ex- 
plain why the equation of state of our DBI dark energy 
approaches —1 in the end. Therefore, our model should 
be treated as a phenomenological model which requires 
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